Measurement of Low-Energy Cosmic-Ray Antiprotons at Solar Minimum 
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The absolute fluxes of the cosmic-ray antiprotons at solar minimum are measured in the energy 
range 0.18 to 1.4 GeV, based on 43 events unambiguously detected in BESS '95 data. The resultant 
energy spectrum appears to be flat below 1 GeV, compatible with a possible admixture of primary 
antiproton component with a soft energy spectrum, while the possibility of secondary antiprotons 
alone explaining the data cannot be excluded with the present accuracy. Further improvement of 
statistical accuracy and extension of the energy range are planned in future BESS flights. 



PACS numbers: 98.70.Sa, 95.85.Ry 

The strict upper limit jy on the cosmic-ray anti- 
helium/helium ratio provides most direct evidence for 
Galaxy and nearby part of Universe being composed 
solely of particles. If so, antiproton (p) in cosmic-ray 
must have been produced in pairs with protons (p's) by 
some elementary-particle processes in Galaxy or in Uni- 
verse. One such process which should certainly exist is 
the interaction of high energy cosmic rays with interstel- 
lar gas. The energy spectrum of p's from this "secondary" 
process is expected to show a characteristic peak at 2 ~ 3 
GeV and sharp decreases of the flux below and above 
the peak. One can also conceive more novel or exotic 
elementary-particle processes such as the annihilation of 
neutralino dark matter or the evaporation of primordial 
black hole Q . The p's from these "primary" sources are 
expected to show very soft energy spectra ||, peaking 
toward the lower energies, and would exhibit large solar 
modulations [Q. 

Although general characteristics of cosmic-ray p data 
accumulated over last 20 years seem compatible with the 
secondary p, poor statistics of the data do not allow 
further investigation. A detector with much larger ac- 
ceptance is needed both to determine the secondary p 
flux and to perform a sensitive search for low-energy pri- 
mary p component. We report here (for detail see ||) 
on such an attempt with the BESS spectrometer, which 
was designed || and constructed for this purpose. By 
using BESS '93 data, the first mass-identified detection 
of cosmic-ray p's was reported QJ^] based on four events 
detected in the low-energy range of 0.3 to 0.5 GeV. This 
measurement (BESS '95) extends the range to 1.4 GeV, 
and was conducted at a period close to solar minimum. 

The BESS detector J§ is shown in Figjl} The thin 
superconducting coil |M (4 g/cm 2 thick including the 
cryostat) produces a uniform axial magnetic field of 1 



Tesla. The r</>-tracking in the central region is performed 
by fitting up to 28 hit-points, each with 200 /im reso- 
lution, in the JET- and IDC-drift chambers, resulting 
in a magnetic-rigidity |ll| resolution of 0.5 % at 1 GV/c. 
Tracking in the z-coordinate is done to an accuracy of 300 
/im by fitting points in IDC measured with vernier pads 
and points in the JET chamber measured using charge- 
division. The continuous and redundant 3-dimensional 
tracking with the drift chambers, all equipped with multi- 
hit capacity, enables us to recognize multi-track events 
and tracks having interactions and scatterings, thus min- 
imizing the background. The dE/dx of the particle in 
the JET chamber is obtained as a truncated mean of the 
integrated charges of the hit-pulses which compose the 
track. The TOF scintillator hodoscopes, newly built for 
'95 flight, measure the time-of-flight of particles with a 
resolution of 110 ps as compared to previous 280 ps. 
The first-level trigger is provided by a coincidence be- 
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FIG. 1. Cross-sectional view of the BESS '95 detector 
with one of the p events. 
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tween the top and the bottom scintillators, with the 
threshold set at 1/3 of the pulse height from vertically 
incident minimum ionizing particles. The second-level 
trigger, which utilizes the IDC and ODC hit-patterns, 
first rejects the null- and multi-track events (pattern- 
selection) and makes a rough rigidity-determination to 
select negatively-charged particles (rigidity-selection) . In 
addition, one of every 90 first-level triggers is recorded, 
irrespective of the second-level trigger condition, to build 
a sample of unbiased triggers from which the efficiencies 
are determined. The '95 BESS flight was carried out on 
July 25, from Lynn Lake, Canada. The scientific data 
were taken at an altitude of 36 km (a residual atmo- 
sphere of 5 g/cm 2 ) with a rigidity cutoff ranging 0.3 to 
0.5 GV/c. During the live time of 2.72 x 10 4 sec, 4.6 mil- 
lion events were recorded on magnetic tapes, out of 60 
million cosmic rays which passed through the detector. 

The following off-line selections are applied equally 
for positive- and negative-rigidity events, (i) Only one 
counter is hit in each layer of the TOF hodoscopes; (ii) 
Only one track, which is fully contained in the fiducial 
region, should be found in the JET chamber. The follow- 
ing cuts are then applied to ensure the qualities of the 
track and the timing measurements. (1) The fitted rib- 
and z-track, respectively, should contain at least 10 and 
5 hits in the JET chamber, and at least one hit in each of 
the upper two and lower two IDC layers; (2) The reduced 
X 2 of both the fitted rcj)- and z-track have to be less than 
5; (3) The extrapolated track should cross the fiducial 
region of TOF scintillators (\z\ < 46.5 cm); (4) The z- 
position determined by the left-right time difference mea- 
sured by the PMTs should match the z-impact point of 
the extrapolated track at the TOF counter within 3.8 cm; 
(5) The ratio of the signal amplitude of the left and right 
PMTs should be consistent with the ^-impact point of 
the extrapolated track. The combined efficiency of these 
off-line selections is found to be 81 % at 0.7 GV/c and 
78 % at 2 GV/c for the protons. 

We then require that the p's and the protons must have 
proper dE/dx in the top and bottom scintillators as well 
as in the JET chamber (shown in FigJ|). These dE/dx 
selections reject all of the particles with charge greater 
than 1 and most of the low-energy e/jj,, while keeping 90 
% efficiency for protons. The /3 _1 vs. rigidity(i?) plot 
for the remaining events is shown in Fig.^. A clear band 
of p's is visible at the position exactly symmetric to the 
proton band. Each of the p candidates was closely inves- 
tigated on the event display for correctness of the track 
fitting etc., and no particular problems were found. The 
probability of positive-rigidity particles or albedo (up- 
going) particles faking any of the p candidates is negligi- 
bly small. Also it is apparent from the cleanness of the 
p-band that there is no background, particularly in the 
high /3 _1 region. We limit further analysis to the region 
> 1.095 in order to have a negligible contamination 
(less than 0.5 events) due to the tail of the e~ / ' \i~ distri- 



bution. The rigidity region is limited to \R\ < 2.1 GV/c, 
since the efficiency of the above cut becomes lower 
than 50 % outside of this |i?|-range. As seen in Fig.| all 
of the 43 negatively-charged events observed in this /3 _1 - 
R region reside inside the 97 % inclusion band expected 
for the p's. Furthermore, the masses determined for these 
events show the correct central value, and the distribu- 
tion is similar to the one observed for the protons. All 
these background-studies and the consistency-checks al- 
low us to conclude that the 43 events are real p's incident 
on the BESS spectrometer. To check against the possibil- 
ity that any of the p's were re-entrant albedo, we numer- 
ically traced back the trajectories of all p's through the 
Earth's geomagnetic field and confirmed that all events 
came from outside of geomagnetic sphere. 

Based on these 43 events, we obtain the p energy spec- 
trum at the top of the atmosphere (TOA) in the follow- 
ing way: The TOA energy of each event is calculated by 
tracing back the particle through the detector material 
and air by using Bethe-Heitler routine of GEANT 3.21. 
The corrections are usually small, 30 ± 3 MeV for a 1 
GeV event. For the particular lowest-energy event which 
has 157 MeV at the center of the detector, the energy 
correction amounts to 66 MeV with estimated accuracy 
of ± 2 MeV. Among the factors necessary to obtain the 
flux, the geometrical acceptance (G) can be calculated 
reliably both analytically and by Monte Carlo methods 
due to the simple geometry and uniform magnetic field. 
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FIG. 2. Scatter plot of dE/dx vs. rigidity (a) for the top 
scintillator and (b) for the JET chamber. The dashed lines 
indicate the dE/dx band for the proton/antiproton selection. 
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FIG. 3. The identification of p events. The solid lines 
define the /3~ 1 -R region and the p mass band used for the 
spectrum measurement. 

The efficiencies of the pattern-selection (e pa t) and of the 
rigidity-selection (e r ; g ) for the second-level trigger as well 
as the off-line selection efficiency for the protons are di- 
rectly determined by using the unbiased trigger sample. 
The selection efficiency for the p's (e se i) is obtained by 
multiplying the proton selection efficiency by the ratio 
(ranging from 1.30 at 0.3 GeV to 1.12 at 1.4 GeV) of 
Monte Carlo simulated selection efficiencies for the p's 
to protons. In this way we minimize the effect of inac- 
curacy of the simulation. The GEANt/gheisha code is 
used for the simulation by incorporating the detailed ma- 
terial distribution and realistic detector performance in 
order to also evaluate the interaction loss in the instru- 
ment (edct)- Care is taken to use the correct p- nuclei 
inelastic (including the annihilation) cross sections, since 
original GHEISHA uses cross sections much higher than re- 
cent measurements, especially below 1 GeV. The survival 
probability (e a i r ) of the p's through the air is estimated 
by the Monte Carlo simulation. We also subtract the ex- 



pected number of atmospheric p's, produced by the colli- 
sions of cosmic rays in the air. Among recent calculations 
]l2| |l4| of the atmospheric p's, which agree within ± 20 
% relative accuracy, we take the one jl2| which utilizes 
the most detailed nuclear model and the 3-dimensional 
Monte Carlo simulation. 

Table |l] summarizes all relevant quantities and resul- 
tant p fluxes in five energy bins. The first and the second 
errors, respectively, represent the statistical and system- 
atic errors. The statistical errors are calculated properly 
p5| for the Poisson distributions. The dominant system- 
atic error at low-energy bins is the uncertainty in the 
interaction losses, to which we attribute ± 40 % rela- 
tive error. At high energy bins, the uncertainty in the 
atmospheric p calculations, to which we attribute ± 30 
% relative error, becomes an important systematic error. 
In all bins, the systematic errors are overwhelmed by the 
statistical errors. Table | also contains the p/p ratios, 
obtained by counting the number of protons which sur- 
vive the same selections as for the p's. The interaction 
losses of the protons in the air and in the instrument are 
corrected by the same Monte Carlo simulations, and the 
atmospheric protons are subtracted based on Ref. fuj ]. 

Our data on the p flux are shown in Figji] together 
with previous measurements and various theoretical cal- 
culations. Our data are consistent with recent measure- 
ments (including BESS '93 result), all of which have much 
larger statistical errors. The energy spectrum we have 
obtained appears to be flat below 1 GeV within present 
statistical accuracy, and does not exhibit the steep de- 



cline which seems generic to calculations 19-21 for the 
secondary p's. This might be due to a statistical fluctua- 
tion, might indicate that the propagation models need to 
be modified, or might perhaps suggest a contribution of 
primary p component with energy spectrum much softer 
than the secondary p. We note [£3| that an admixture 
of p's [|| from evaporating primordial black holes (PBH) 
0, added to a standard leaky box (SLB) calculation ||l]] 
of the secondary p's, could provide a fit (thick dashed 
curve in Fig^) to the data. Primary p's component from 
the annihilation of ncutralino dark matter could also ex- 
plain |^,^3| the observed spectrum, if we assume an astro- 
physical enhancement factor of 15 or more, which might 
be possible jl| due to the dumpiness of the dark mat- 
ter. However, given the present large statistical errors 
and the ambiguities in the flux of the secondary p's, we 



TABLE I. Summary of antiproton flux measurement. 



TOA Kinetic energy 


No. of 


No. of G 


Epat trig Cscl Edct Eair TOA p flux 


TOA p/p ratio 


(GeV) 


detected p 


atm. p (m 2 sr) 


(m-W^GeV" 1 ) 

n-rn n nc n co r\An nan 1 ne+1.41 +U.3U , „-2 


n -7D+U.81 +0.18 m-5 



0.5 
0.7 



0.3 
0.5 
0.7 - 1.0 
1.0 - 1.4 



7 
7 
11 
15 



0.78 
1.4 
2.8 
3.5 



0.31 
0.31 
0.32 
0.32 



0.73 
0.73 
0.72 
0.72 



0.92 
0.87 
0.80 
0.72 



0.74 
0.74 
0.73 
0.56 



u.^a u.oi r.oo_ 79 _ 30 x iu 

0.62 0.88 1.36±g;|f ±g;li x 10~ 2 

0.65 0.89 l-22+°;®5l°;2o x 10 -2 

0.67 0.89 1.25to:5ota22 x 10 ~ 2 

o.68 o.9o i.85±g:^f ±S:li x 10-2 



'-0.45 -0.18 
n 74 +0.45+0.12 ln -5 

u - '^-0.31 -0.12 x iu 

+0.52+0.12 ln -5 

U - ' '-0.35 -0.12 X 1U 

+0. 54+0.17 ln -5 

1 - ul -0.40 -0.17 x lu 

1 qq+0 85 +0.32 ln -5 

±.»»_ 66 _ .32 A lu 
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cannot exclude the possibility of the secondary p's alone 
explaining the data, especially if the actual secondary p 
flux is among the largest predicted. To clarify the situ- 
ation, one needs; a) to improve the statistical accuracy, 
b) to detect the predicted peak of the secondary p's by 
extending the energy range of the identification and to 
accurately measure its absolute flux, and c) to measure 
the change of the energy spectrum with solar activity. 
Concerning the point c), it was shown B that a clear 
signature of a primary p component can be obtained by 
observing an enhancement of low-energy p's at the solar 
minimum and a sharp decrease afterward. Independent 
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FIG. 4. Comparison of the BESS '95 p fluxes with pre- 
vious data and various calculations. The upper (1) and the 
lower (2) thin lines represent respectively bounds in standard 
leaky box calculation by Simon et al. [19]. The middle thin 
line (3) is the median flux calculated in the standard leaky box 
model by Gaisser and Schaefer [20], who attribute a factor 2 
ambiguities in both directions. The upper (4) and the lower 
(5) thick lines are the calculations by Mitsui et al. [21] in 
the standard leaky box model and in the diffusive reaccelera- 
tion model (DR) [22], respectively. The thin (6) and thick (7) 
dashed lines [23] are, respectively, the spectrum from evap- 
orating primordial black holes (with local explosion rate of 
6 x 10 -3 pc _3 yr _1 ) and its sum with the standard leaky box 
model calculation by Mitsui et al. The thin (8) and thick 
(9) dotted lines, respectively, are the spectrum from the an- 
nihilations of neutralino dark matter [4] and its sum with the 
diffusive reacceleration calculation by Mitsui et al. All theo- 
retical fluxes are modulated [24] with (f>F = 530 MV, which 
corresponds to the BESS '95 flight. 



to the search for the primary p component, the precise 
measurement of the secondary p spectrum itself will be of 
crucial importance to determine the propagation mecha- 
nism of cosmic rays in the Galaxy. We intend to accom- 
plish these purposes in future BESS flights. 
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